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1 Introduction 

The shell model is a fundamental many-body approach to the study of atomic 
nuclei [||]. It explains the magic numbers as shell closures and the energy 
spectra of odd-mass nuclei near closed shells as that of the odd nucleon 
in a potential well defined by the closed shell nucleons. The remarkable 
advances in computer power and the use of complex algorithms have allowed 
for systematic studies of most sd- and fp-sheW nuclei [§]. However, in heavy 
nuclei is not possible to solve the shell-model problem exactly. Although the 
fermionic character of the nucleons restricts their allowed degrees of freedom, 
the number of accessible states of a system still grows combinatorially with 
the number of valence nucleons. For this reason truncation schemes must be 
introduced. 

In light deformed nuclei the dominance of the quadrupole-quadrupole in- 
teraction led to an introduction of the SU(3) shell model 0. Within the 
SU(3) algebraic framework, large shell-model spaces can be truncated in a 
very natural way. While in general realistic interactions mix irreducible repre- 
sentations (irreps) of SU(3), the ground state wavefunction of well-deformed 
light nuclei are typically dominated by a few SU(3) irreps @, ^ |^. The 
strong spin-orbit interaction renders the SU(3) truncation scheme useless in 
the heavier nuclei of the fp-shell, while at the same time for rare earth and 
actinide species pseudo-spin emerges as a good symmetry and with it the 
pseudo-SU(3) model 0, |, i. 

Pseudo-spin can be recognized from the experimental fact that single- 
particle orbitals with j = I — 1/2 and j = {I — 2) + 1/2 in the shell r) lie 
very close in energy and can therefore be labeled as pseudo spin doublets 
with quantum numbers ] = j, rj = r] — 1, and / = / — 1. The origin of 
this symmetry has been traced back to the relativistic mean field equations 

In this work the energy spectra and B(E2) transition strengths of ^^^Tb 
are calculated from a microscopic perspective using the pseudo SU(3) shell 
model. Normal parity bands are described in a many-body basis built with 
active nucleons occupying normal parity levels. Polarization effects due to 
valence nucleons in intruder orbits are taken into account through the use of 
effective charges. 

Basis states are built from SU(3) irreps obtained by taking the direct 
product of proton and neutron representations. In |l^, ^ it is shown that for 
a description of the low-energy spectrum of deformed nuclei, the Hilbert 



space should be truncated according to contributions of the quadrupole- 
quadrupole interaction and the single-particle one-body Hamiltonian. While 
pairing plays a very important role in determining the energy spectrum, but 
does not strongly modify the wave functions of deformed nuclei. 

As can be seen in |1T4| , ^, by using an schematic Hamiltonian parametrized 
according to systematics [|16], ^7^, it is possible to describe the low-lying en- 



ergy spectrum of even- and odd-mass heavy deformed nuclei. While the 
application of the model to other deformed rare-earth and actinide nuclei is 
in order, in the present contribution we will address some specific questions 
about the pseudo SU(3) model. In particular, we will discuss the importance 
of various terms in the Hamiltonian, ways in which their strengths can be 
deduced from systematic, and the effect each term has on the energy spec- 
tra. Taking ^^^Tb as an example, we will examine the wave functions of each 
rotational band, paying particular attention to the B{E2] J ^ [J — 1)} and 
B{E2] J ^ (J — 2)} transition strengths between states in the same bands 
and between states belonging to different bands. 

In section 2 the pseudo SU(3) classification scheme is presented. The 
pseudo SU(3) Hamiltonian and its parametrization is discussed in Section 3. 
The effect which each term in the Hamiltonian has on the energy spectra is 
analyzed in Section 4. Section 5 contains the analysis of the wave functions 
associated with different rotational bands and the associated B(E2) values. 
Conclusions are discussed in Section 6. 

2 The pseudo SU(3) basis 

The first step in any application of the pseudo SU(3) model is to build the 
many-body basis. To do this it is necessary to know how many valence 
nucleons occupy the normal parity orbitals. We will show how this is done 
using ^^^Tb as an example. It has 65 protons and 94 neutrons, and of these, 
15 protons and 12 neutrons are in the last unfilled (open) shells. Assuming a 
deformation (3 ~ 0.25, the deformed Nilsson single-particle levels of the active 



shells are filled from below |18, IS]. Nine protons are distributed in the lg7/2 
and 2(^5/2 orbitals of the rj = A shell, and the remaining six occupy the l/iii/2 
intruder orbital. Eight neutrons occupy the 2/7/2 and 1/19/2 orbitals of the 
f] = 5 shell and four are in 1^13/2 orbital. The relevant occupation numbers 
nT^^riy can be summarized as 



n^ = 9, < = 6, < = 8, n^ = A (1) 

where N refers to normal parity states and A to the abnormal parity (also 
called unique or intruder) states. The deformed Nilsson mean field is only 
employed to define the number of nucleons in normal and unique parity 
orbitals. The use of the pseudo SU(3) basis to describe the normal parity 
sector implies that these nucleons can occupy all normal parity orbitals, not 
only those with the lower single particle energies. 

As it has been the case for all pseudo SU(3) studies up to now, we will 
freeze the nucleons in abnormal parity orbital and describe the dynamics 
using only nucleons in normal parity states. While it has been shown that 
this is a reasonable, it is nonetheless a strong assumption. This choice is 
further reflected through the use of effective charges to describe quadrupole 
electromagnetic transitions which are larger than those usually employed 
in typical shell-model calculations for light nuclei. A more sophisticated 
treatment of the problem, with nucleons in intruder orbitals described in the 
same footing using SU(3) irreps is under development 0. 

The many-particle states of Ua active nucleons in a given normal parity 
shell rjai a = V oi Ti can be classified by the following chains of groups: 



{!""} {fa] {fa] la{K,^^a) Sa Ka 

U{n^) D U{n^/2) X f/(2) D SU{3) x SU{2) D 

^a ^Q 

50(3) X SU{2) D SUj{2), (2) 

where above each group the quantum numbers that characterize its irreps 
are given and ja and Ka are multiplicity labels of the indicated reductions. 

Any state | JjM), where J is the total angular momentum, M its projec- 
tion and i an integer index which enumerates the states with the same J, M 
starting from the one with the lowest energy, is built as a linear combination 

\J,M) = J2C'p^''\f3JM) (3) 

of the strong coupled proton-neutron states 

\(3JM) = \{f\]{X^fi^)S^,{U]{X,fi,)S,;p{Xfi)KL,S JM) 



= Y. (LMl^SMsIJM) Y. iS^Ms.,S,Msu\SMs) 

MlMs MsnMs^ 

/vTT Ai.j/ Ijty l^U IVlTZ IvliJ 

In the above expression (— ; — |— ) and (— , — |— ) are the SU(3) and SU(2) 
Clebsch Gordan coefficients, respectively. We are considering only configu- 
rations with the highest spatial symmetry [^, |20|. For ^^^Tb the active shells 



in the pseudo SU(3) space are f;^ = 3 and r]y = A with degeneracies ^2^ = 20 
and Vty = 30, respectively. In the large groups U{10) and f/(15), the spa- 
tially most symmetric irreps for 9 protons and 8 neutrons are, respectively, 
{U} = {2^1} and {/^} = {2^}. It implies that S^ = 1/2 and S^ = 0. In 
other words, we are only taking into account configurations with pseudo spin 
zero for an even number of nucleons and 1/2 for an odd number of nucleons. 

The above considerations rely strongly on the goodness of pseudo-spin 
symmetry, which manifests itself in the near degeneracy of the pseudo spin- 
orbit partners. When the real SU(3) model is used to describe deformed 
nuclei in the p/-shell, the spin mixing is very important and Eq. (1) must be 
modified accordingly |^, ^ . 

What makes the pseudo SU(3) model a powerful theory is that it al- 
lows one to envoke a relatively simple and physically motivated basis trunca- 
tion scheme. Extended shell-model calculations in the pf- and sc^^^shell have 
shown that in the description of deformed nuclei the Hilbert space can be 
truncated to only those states that are relevant when both the quadrupole- 
quadrupole force and the single-particle Hamiltonian are taken into account 
]TB| . While pairing is fundamental to obtaining the correct moment of inertia 
of the rotational bands, it has a relatively small effect on the overall wave 
functions \T^\. An analysis of the SU(3) content of wave functions obtained 



in large shell- model diagonalizations |21|], as well as the excellent description 
of ground and excited bands in heavy deformed even- n, |22[ and odd-mass 



| T5| nuclei strongly support this statement. 

The quadrupole-quadrupole interaction can be expressed in terms of the 
second order SU(3) Casimir operator C2, 

Q.g = 4C2-3L(L + l). (5) 

The eigenvalue of C2 for a given of SU(3) irrep (A, /i) is given by 

(C2) = (A2 + /i2 + A/i + 3A + 3/i). (6) 



The larger the expectation value of C2, the greater the binding of that SU(3) 
irrep by a pure Q ■ Q interaction. We build the basis selecting the proton 
and neutron irreps with the largest (C2). 

For ^^^Tb, the SU(3) representations were selected as follows. Including all 
the possible spins there are 97 irreps in the proton space and 285 irreps in the 
neutron space. Proton and neutron irreps which belongs to the irreps {/tt} = 
{2^1} and {fi^} = {2^} of the large groups f/(10) and f/(15), respectively, 
are found and ordered according to their C2 value. In Table 1 the seven 
representations with the largest C2 in ^^^Tb are shown. 

Table 1 

Taking the direct product of these protons and neutrons irreps results in 
many strong coupled SU(3) irreps. From these, the 15 with the largest C2 
were chosen. They are shown in Table 2. These 15 proton-neutron irreps 
define the Hilbert space of the model. They are a small subset of all the pos- 
sible irreps, and involve only 4 proton and 3 neutron irreps. In this strongly 
truncated space it is possible to describe the low energy spectra of even-even 
nuclei |I^ p2| and normal parity bands in odd- mass nuclei |T^. The va- 



lidity of the pseudo SU(3) symmetry is the rationale behind this successful 
truncation scheme. 
Table 2. 

3 The pseudo SU(3) Hamiltonian 

The Hamiltonian contains spherical Nilsson single-particle terms for protons 
(Hsp^Tr) and neutrons (Hsp^u), the quadrupole-quadrupole (Q • Q) and pairing 
interactions {Hpair,TT and Hpair,u), as well as three 'rotor-like' terms which are 
diagonal in the SU(3) basis. 

H = 22 {Hsp,a — Ga Hpair,a} — 7^ X Q ' Q (7) 

a=n,u ^ 

+ aKJ + b J^ + Aasym Cl- 

This Hamiltonian can be separated into two parts: the first row includes 
Nilsson single-particle energies and the pairing and quadrupole-quadrupole 
interactions (Q is the quadrupole operator in the pseudo SU(3) space, see 
below). They are the basic components of any realistic Hamiltonian [|1^, [1^ 



and have been widely studied in the nuclear physics literature, allowing their 



respective strengths to be fixed by systematics |TB, 0. In the second row 
there are three rotor terms used to fine tune the moment of inertia and the 
position of the different K bands. The SU(3) mixing is due to the single- 
particle and pairing terms. 

The three 'rotor-like' terms have been studied in detail in previous papers 



where the pseudo SU(3) symmetry was used as a dynamical symmetry |lg, 
|T9| . In the present work, a and h are the only two parameters used to fit the 
spectra. 

The term proportional to K^ breaks the SU(3) degeneracy of the different 
K bands |^. It has the form 

^2 ^ A1A2J + XzX^ + XI .. 

/A3 + A1A2 
where J^, X| and X| are the three rotational scalars formed with products of 



J and Q and the Aj coefficients are functions of A and /i pO| , p!8| , p!9|] , as shown 
in Appendix A. The reduced matrix elements for X3 and X^ are evaluated 
using Racah and SU(3) coupling coefficients [Q |2^. For ^^^Tb a = 0.0198 
was found to provide the best fit. 

The term proportional to J^ is used to fine tune the moment of inertia. 
It represents a small correction to the quadrupole-quadrupole term, which 
contributes to the rotor spectra with strength 3/2% (see Eq. (^)). For ^^^Tb 
we used b = —0.0031, which introduce a change of about 15% in the rotational 
spectra. 

The asymmetry term distinguishes SU(3) irreps with both A and /i even 



from the others ||2^, having no interaction strength in the first case and a 
positive one in the the second. In this way the contribution of irreps with both 
A and /i even is slightly enhanced because they belong to different symmetry 
types of the intrinsic Vierergruppe D2 |2^. The asymmetry coefficient has 
a value Aasym = 0.0008, fixed according with ||2^. The same value was 



employed for the three A=159 nuclei studied in |T3 
The single-particle Nilsson Hamiltonian is 



Hsp = hujoi^r] + ^) - KhuJo{2l ■ s + fil }, (9) 



with parameters |Ti 

hcuo = 4:lA'^/^[MeV], k^ = 0.0637, k^ = 0.0637, (10) 



/i^ = 0.60, /i^ = 0.42, 
The pairing interaction is 

Vp = —-G'^Ojalaj/ajr (11) 

jd' 

where J denotes the time reversed partner of the single-particle state j and G 
is the strength of the pairing force. Its second quantized expression in term 
of SU(3) tensors is reviewed in Appendix B. For the pairing coefficients Gtt.i/, 
we used [jl6|, |17| 



G^ = ^- = 0.132, G^ = L = 0.106. (12) 

In the SU(3) model the collective quadrupole operator, defined by Q^ = 



167r/5X]j^j^^2/x(^i)/&^5 is symmetrized in order to correspond to one of the 
SU(3) generators. It is called the 'algebraic' quadrupole operator Q'^ = 

4:Ti/5J2i[rfY2^{fi)/b'^ + b'^pfY2f_i{pi)] ||18|, [l9l. Within a major oscillator shell. 



the matrix elements of Q'^ and Q'^ are identical. When transformed to the 
pseudo SU(3) basis, it maps to a linear combination of SU(3) tensors which 
is dominated by the quadrupole operator Q, which is the generator of the 
pseudo SU(3) algebra [|1^. It is this quadrupole operator in the pseudo 



SU(3) space the one included in the quadrupole-quadrupole interaction in 
Hamiltonian (^. In first order, the relationship between both quadrupole 
operators can be written as 

Q", = ^Qm (13) 

which holds for protons and neutrons separately. 
The coefficient x of the operator Q ■ Q is 

35 
X = ^^ = 0.00753. (14) 



It is consistent with the parametrization discussed in |T^, provided one keep 
in mind that in this reference the quadrupole operator is just r^y2^; and 
that, as mentioned above, when operating in the pseudo SU(3) space the 

interaction must be, in first order, a factor (^^) stronger than the similar 



one in the normal space. This imphes that the constant x i^i 1131 ^^ ^ factor 



167r/5 
[('7+1)/^] = 



6.5 larger than our x- 



The electric quadrupole operator is expressed aspO 



C^/i G-jT^TT ~r Gy^y ^'^ 6-71- ^:7^ ~r ^y ^U} \ ^/ 

with effective charges e,r = 2.3, e^ = 1.3. These values are very similar to 
those used in the pseudo SU(3) description of even-even nuclei [|^, |14|. They 



are larger than those used in standard calculations of B(E2) strengths |T^ 
due to the passive role assigned to the nucleons in unique parity orbitals, 
whose contribution to the quadrupole moments is parametrized in this way. 

4 The energy spectra 



Fig. 1 shows the yrast and excited bands in ^^^Tb. Experimental data ||2^ 
are plotted on the left hand side, while those obtained using the Hilbert space 
and the Hamiltonian parameters discussed in the previous sections are shown 
in the right hand side. The agreement between both is excellent, as is the 
case for the nuclei ^^^Eu and ^^^Dy, whose energy spectra was studied in |T5[ . 
From the four bands reported in the literature in ^^^Tb, three of them (the 
yrast, 5/2-^ , and 3/22 bands) have a difference between the experimental and 
predicted levels of less than 50 KeV. The 1/2^ is slightly high in energy, and 
the model predicts an exaggerated staggering, with origin discussed below. 

Figure 1 

The whole energy spectra is built up by the interplay between the single- 
particle and quadrupole-quadrupole terms in the Hamiltonian |TB, ^. These 



two terms define the relative ordering between the different bands, as well as 
the main components of the wavefunction. As expected, the use of realistic 
single-particle energies plays a key role in the appropriate description of odd- 
mass nuclei. 

To make this point clear, in Fig. 2 the theoretical energy spectra calcu- 
lated with Hamiltonian (|^) without single-particle energies are presented on 
the right of each column and compare it with the corresponding experimental 
energies on the left of each column |2^. It is clear that the ordering is shifted. 



The ground state is now predicted to have J = |, the first excited band starts 
with J = |. Only the third band with J = \ reamains ordered correctly. 
The wave functions of all the states in all these three bands are dominated 



(more than 95%) by the leading SU(3) irrep (28,8), (10,4)^(18,4)^. We will 
return to this point when discussing the wave functions in Section 5. 

Figure 2 

From Fig. 2 it is also clear that, in absence of single-particle energies, the 
quadrupole-quadrupole interaction rules over all the others, even in the pres- 
ence of realistic pairing strengths. It determines the single SU(3) irrep that 
dominates the low energy spectra and the rotor pattern each band exhibits. 

In order to discuss the effect each of the remaining term in Hamiltonian 
d^) has on the energy spectra, three rotational bands are plotted in Figure 
3: the yrast band in insert (a), the 5/2^ band in insert (b) and the 1/2^ 
band in insert (c). In the first column from the left the experimental en- 
ergies are plotted for each band, with the angular momentum and parity 
of each state written on the left. The next column presents the theoretical 
results, the same ones shown in Fig. 1. The next three columns show the 
theoretical spectra obtained without the pairing interaction (third column, 
G^ = Gj, = 0), without pairing between protons (fourth column, G^ = 0) 
and without pairing between neutrons (fifth column, G^, = 0). The remain- 
ing four columns depict the behavior of the spectra when the rotor terms are 
turned off. The sixth column, labeled K"^, shows the spectra without this 
term (a = 0), the seventh column without the J^ term (6 = 0), the eighth 
without the asymmetry term {Aasym = 0) and the ninth presents the spectra 
with all the three rotor terms turned off (a = 6 = Aasym = 0). The last 
column on the right presents the experimental energies again, to help with 
the comparison. 

Figure 3 

The theoretical description of the yrast band is good. The effect of the 
pairing interaction is clearly seen in the third column; namely, it expands the 
energy spectra. This is a rather remarkable result that seems to comes about 
because of the highly truncated nature of the basis. The pairing interaction 
strongly mixes the SU(3) irreps [^ ^ ^ ^] and plays an important role 



in determining the moment of inertia of this deformed nuclei. At the same 
time its effects on the wavefunction are minor, as has been discovered in 
indenpendent shell-model calculations [^ . Comparing the third, fourth and 
fifth columns, it is clear it is the pairing interaction between neutrons that is 
most important. The last uncoupled proton in ^^^Tb seems to nearly ignores 
the effect of pairing. The K^ term has negligible effect on the yrast band. 
The J^ term helps to fine tuning the moment of inertia to the experimental 
number. The asymmetry term pushes the J = y state down in energy and 

10 



closer to the experimental value. Although a small contribution it is clearly 
important and provides justification for including this term included in the 
Hamiltonian. 

The 5/2j^ band shown in insert (b) shows a similar behavior. The first 
three states reported are in close agreement with the experimental data, 
and theory predicts more states that belongs to this band. In the absence 
of pairing, the energy spectra is compressed which is due mainly to the 
protons. The K^ term is clearly relevant for this band: when it is missing 
the band head moves down about 100 KeV in energy. The J^ term modifies 
the moment of inertia slightly. 

The 1/2-,^ band is shown in insert (c). Its band head is predicted at an en- 
ergy 150 KeV larger than the measured one. The theory predicts three nearly 
degenerated pairs of levels: (|, |), (|, |) and (y' ^) while their experimental 
counterparts are less closely packed. As can be seen in the eighth column, 
when Aasym = the level spacing is closer to the reported one, but the band- 
head energy it too high. It is noticeable than turning off the K^ term the 
energies move upwards, while in the 5/2^ band they move downward. 

The contribution of each term in the Hamiltonian can be summarized as 
follows: 

• Quadrupole-quadrupole: This interaction is tied to the quadrupole 
deformation, and since the pioneering work of Elliott [Q] has been known 
to play a crucial role in the dynamics of deformed nuclei. In the pseudo 
SU(3) model its dominance can also be used to determine the trunca- 
tion of the Hilbert space. 

• Spherical single-particle energies: These form the basis of the shell 
model, and serve define the low-energy spectra of odd-mass nuclei. And 
this is precisely the role they play in the pseudo SU(3) model: they de- 
termine to a large extent relative ordering of the different bands. To- 
gether with the quadrupole-quadrupole interaction, the single-particle 
energies define the gross features of the energy spectra and the mixing 
of different SU(3) irreps in the wave functions. 

• Pairing: This interaction expands the whole energy spectra, almost as 
though it is a multiplicative constant in front of Hamiltonian (^. It also 
shifts the energy of each band head and alters the moment of inertia of 
each band. The success of previous investigations that used the pseudo 



SU(3) symmetry as a dynamical symmetry [18, 19| without the mixing 
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of SU(3) representations, is partially a result of this effect: a larger 
quadrupole-quadrupole and rotor strength mimics the effect of pairing 
on the spectra. The pairing interaction acting in the subspace with an 
odd number of nucleons has a nearly negligible effect on the low-lying 
bands. The subspace with even number of nucleons is responsible for 
most of the pairing effects. 

K^: Since it has been shown that this particular combination of 2-, 3- 
and 4-body terms corresponds to the square of the third component of 



the angular momentum in the intrinsic frame p5[, it can be used to 
adjust band-head energies. For the 5/2-^ band, when a = the band- 
head energy is 277 keV, while in the presence of this term it moves up 
to 356 keV, which is close to the experimental value of 348 keV [^. 

• J^: This diagonal term provides small corrections to the moment of 
inertia. The negative b value used in the present study makes serves to 
increase the moment of inertia, compressing the corresponding spectra. 

• Asymmetry: This term enhances the contribution of the SU(3) irreps 
with both A and fi even. It has an important effect on states with large 
angular momentum. For example, in the absence of this term the state 
Y of the yrast band is displaced to higher energies due to mixing with 
other SU(3) irreps with A or /i odd. 

• Rotor terms: The effect of suppressing simultaneously the three rotor 
terms, i.e. taking a = b = Aasym = 0, is shown in the ninth column of 
Fig. 3. The inclusion of Nilsson single-particle energies and the pair- 
ing and quadrupole-quadrupole interactions suffices to provide a very 
reasonable energy spectra, with all known bands in their correct order 
and the overall rotor features reproduced. The rotor terms provide the 
fine tuning, with energies being adjusted by no more than 15 %. The 
predictive power of the pseudo SU(3) model strongly relies in this fact. 

5 B(E2) transition strengths and wave func- 
tions 

Up to this point we have centered the discussion on the energetics. But 
the pseudo SU(3) model is far more powerful than this, it can also success- 
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fully describe the electromagnetic transitions. Most of the B(E2) transition 
strengths between states in the first four bands in ^^^Tb are presented in Fig. 
4 for the J — * J — 1 transitions, and in Fig. 5 for the J —>■ J — 2 transitions. 

Figure 4 

Figure 5 

Calculated B(E2) values are given in units of e^6^ x 10^^. They are written 
close to the arrow which graphically show the transition they describe. The 
effective charge employed was discussed in relation with Eq. (pISl). This is the 
only extra parameter introduced, and has a value very close to those used in 
previous studies [^]. To assess the quality of the results, the experimental 



B(E2) values |26| between yrast band states are reported in Table 3. Most of 
the transition strengths are reproduced within the experimental error bars. 
One exception is the transition 17/2^^ — »• 13/2^^ which is underestimated. This 
could be related with the change in the wavefunction of the first y? whose 
mixing with the second ^ state seems to be exaggerated in the model, as 
discussed below. 

Table 3 

Notice that while the intraband B(E2) transition strengths are on the or- 
der of hundreds e^b"^ x 10~^, the interband transitions are much less, typically 
on the order of e^6^ x 10^^ or fractions thereof. This fact supports the identi- 
fication of states belonging to bands, and is consistent with the wavefunction 
analysis. The only measured interband transition is the 1/2^^ -^ 5/2-^ and it 
too is well reproduced by the model, wich predicts a value of 2.96 e^6^ x 10~^. 

It is very instructive to analyze the wave functions of the states belonging 
to the lowest lying energy bands obtained with the Hamiltonian (0). The nine 
most important SU(3) irreps which are relevant to the description of these 
bands are listed in Table 4. 

Table 4 

In Figure 6 the percentage each irrep contributes to each state is plotted 
as a function of the angular momentum of the state for the different bands. 
These were calculated from the wavefunctions, Eq. (0), as simply 100 x 
|^JM«|2^ The symbols listed in the first column of Table 4 identifies the 
various components. All contributions larger than 2 % are plotted, and in 
all cases the states shown add up to at least 95% of the total wavefunction. 

Figure 6 

Each insert, (a) to (d) in Figure 6, gives the main components of one of 
the bands plotted in Figs. 4 and 5, from the band head and up to the state 
with J = ^ or ^. As shown, all have a very regular structure as one moves 



13 



up the bands. While in all the cases there is strong mixing of SU(3) irreps, 
the mixing remains nearly the same for the states with different angular 
moments belonging to the same band. In this sense the mixing is adiabatic 
within each band. This coherence explains the large B(E2) values for intra- 
band transitions. 

In insert (a) the components of each state belonging to the 3/2-^ band 
are shown. As shown, the states from J'^ = 3/2 to 17/2 are about 30 % 
(30,4)[(10,4),x(20,0)J,29%(28,8)[(10,4)^x(18,4)J,14%(30,4)[(ll,2)^x 
(18,4)J, 12 % (31,2)[(11,2)^ X (20,0)J with smaller contributions from the 
(31, 2) and (29, 6) representations that are almost constant for all the states 
in the band. 

States in the 5/2^ band, insert (b), have about 36 % (28, 8) [(10, 4)^ x 
(18,4)J, 22 % (30,4)[(11,2)^ x (18,4)J, and 20 % (30, 4) [(10, 4)^ x (20, 0) J 
which decreases to less then 15% for the J = y state. Other irreps con- 
tribute less than 10% each. The 1/2^ band, insert (c), is the purest of 
the four considered. It has around 60 % (28, 8)[(10,4)^ x (18,4)J, 20% 







(30,4)[(11,2)^ X (18,4)J, and 10 =% (30,4)[(10,4)^ x (20,0)J. It is in- 
teresting to note that this band becomes the ground state band when the 
single-particle energies are not present (see the discussion below Figure 2), 
and in that case all the states in the four bands are built primarily out of 
the (28, 8) [(10, 4)^ x (18, 4) J irrep. 

The 3/22 band, insert (d), is dominated by the irrep (30, 4) [(10, 4)^ x 
(18, 4)^] which get strongly mixed with the (28, 8) [(10, 4)^ x (18, A) J for J = 
^ and f . 

As was mentioned above, the interplay between the single-particle and the 
quadrupole-quadrupole terms in the Hamiltonian defines the SU(3) mixing 
in the wave functions. The four bands discussed here have strong mixing, 
and the ground state band is not dominated by the leading irrep, the one 
which would constitute the ground band in the pure SU(3) symmetry limit. 
The delicate balance between these two interactions, whose strengths are 
taken from known systematics and not used as fitting parameters, defines 
the gross features of the calculated energy spectra which are found to be in 
good agreement with the available experimental information. 

Although the basis is strongly truncated, being built from the 15 SU(3) 
irreps listed in Table 2, not even all of these play an important role. The 
low-lying energy bands discussed above are dominated by the irreps coming 
from the first four rows in Table 2, which are combinations of the two proton 
SU(3) irreps (10,4)^ and (11,2)^ and the two neutron SU(3) irreps (18,4)^ 
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and (20,0)^. 

6 Conclusions 

We have shown that the normal parity bands in the odd-mass heavy deformed 
nuclei with A = 159 can be described quantitatively using the pseudo SU(3) 



model [|T5|. A careful examination of the wave functions and B(E2) transition 
strengths for the four low-lying energy bands in ^^^Tb was made, analyzing 
their structure in terms of their SU(3) components, and relating them with 
their intra- and inter-band transitions, for which the calculated values agree 
closely with the known experimental numbers. 

The most relevant feature of the present application of the pseudo SU(3) 
model is a determination of the primary features of the energy spectra of the 
normal parity rotational bands in ^^^Tb using a Hamiltonian with Nilsson 
single-particle energies and quadrupole-quadrupole and pairing interactions 
with strengths fixed by systematics - strengths of the primary interactions 
were not varied to obtain a "best fit" to the date. A few extra rotor-like 
terms were used to obtain a more precise description of the energies and 
B(E2) values, but this "fine tuning" not affect the spectra in a major way 
and had little influence on the structure of the calculated wavefunctions. 

This work shows the pseudo SU(3) model to be a powerful shell-model 
scheme, one that can be used to describe normal parity bands in deformed 
rare-earth and actinide isotopes by performing a symmetry dictated trun- 
cation of the Hilbert space and using a systematic parametrization of the 
dominant terms in the Hamiltonian. It opens up the possibility of a detailed 
microscopic analysis of other nuclear properties of heavy deformed nuclei, 
both with even and odd numbers of protons and neutrons, like g-factors. Ml 
transitions and beta decays. 
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The residual interaction, K^, is a linear combination of J^, X^ and X4, 
which are rotational scalar operators built from generators of the SU{'i) 
algebra, that is p 



T9[ 



J2 = 


eU? 


^3 = 


Ei,j '^i'^ij'^j 


x,= 


^,,k J^QtjQ^kJk 



(16) 



(17) 

where Jj and Q"- are cartesian forms of the total angular momentum and 
the quadrupole operators, respectively. The K is interpreted to be the third 
component of the total angular momentum along the intrinsic body-fixed 



symmetry axis of the system, which is given by |18, 19 



K^ = (A1A2 J' + A3X3 + X4)/(2A^ + A1A2 



(18) 



with the parameters Aj denoting the eigenvalues of the mass quadrupole 
operator, which are related to the SU{3) labels (A, fj,) through the expressions 



31, 23 



Ai = ^(/i-A), A2 = -^(A + 2/i + 3), A3 = i(2A + /i + 3). (19) 

The last expressions can be obtained by requiring a linear correspondence 
between the invariants of the SU{3) and the semidirect product T5 A 5*0(3) 



groups 1^, ^. Indeed, the expectation value of the operator (|Tq) , with re- 
spect to an orthonormalized basis associated to the chain of groups SU{3) — >■ 
SO {3), corresponds to the eigenvalues of the third component of an intrinsic 
angular momentum when L <^ min{X,^). 

Appendix B: The pairing interaction 

The pairing term (p^Tf ) in the Hamiltonian can be expressed in second 



quantization as p^, E8| 



G 

^P="^ Yl II^»?r,'{(Al,/il)(A2,/i2)Po(Ao,/io)} 



H^alf'^ 



[ttrf' (x> a^/j 



^^ Po>^ofJ-oko=ilo=so=0 



(20) 
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where 



P^^/{(Ai, /ii)(A2, /i2)po(Ao, /io) 



Y: V(2^ + 1)(2^' + 1)((^, 0)/; iv, 0)/| |(Ai, /xi)10) (21) 
w 

((V,0)/';(V,0)/'||(A2,/i2)10)((Ai,/ii)10;(A2,/i2)10||(Ao,/io)10)}. (22) 
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Table Captions 

Table 1: Irreps and C2 values for protons and neutrons in ^^^Tb. Only the 
seven irreps with largest C2 values are listed . 

Table 2: The 15 pseudo SU(3) irreps used in the description of ^^^Tb bands. 

Table 3: Experimental B(E2) transition strengths for ^^^Tb. 

Table 4: Explicit form of the irreps referred to as the components of the 
wavef unctions in Figure 6. 
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{K,fJ'n) 


C2 


{K,fJ^u) 


C2 


(10,4) 


198 


(18,4) 


478 


(7,7) 


189 


(20,0) 


460 


(11,2) 


186 


(16,5) 


424 


(2,11) 


186 


(17,3) 


409 


(8,5) 


168 


(18,1) 


400 


(5,8) 


168 


(13,8) 


400 


(9,3) 


153 


(14,6) 


376 



Table 1 



(A^,/i^) 


{K,fj'u) 






total (A,;u) 


(10,4) 


(18,4) 


(28,8) 


(29,6) 


(30,4) (31,2) (32,0) (26,9) 


(11,2) 


(18,4) 


(29,6) 


(30,4) 


(31,2) 


(10,4) 


(20,0) 


(30,4) 






(11,2) 


(20,0) 


(31,2) 






(7,7) 


(18,4) 


(25,11) 


(26,9) 




(10,4) 


(16,5) 


(26,9) 






(8,5) 


(18,4) 


(26,9) 







Table 2 
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Experimental B(E2) for ^^^Tb 
Jt^{J- 2)i Exp [6^62 X 10- 



7/27- 


-3/2^ 


73.7 ±10.2 


9/27- 


-5/27 


111.5 ±6.1 


11/2+ 


-7/2+ 


147.9 ±3.6 


13/2+ 


-9/2+ 


166.3 ±4.6 


15/27 - 


- 11/27 


158.1 ±10.7 


17/27 - 


-^ 13/27 


154.5 ±10.2 


1/27- 


-5/2+ 


2.56 


4->( 


J -I)) 


Exp [e%^ X 10-2; 


5/27- 


-3/2- 


186.7 ±9.7 


7/27- 


-5/2+ 


117.7 ±25.6 


9/27- 


-7/2+ 


59.9 ±6.6 


11/2+ 


-9/2+ 


57.3 ±5.6 


13/27 - 


- 11/27 


32.2 ±4.1 


15/27 - 


-^ 13/27 


38.4 ±6.6 



Table 3 





labels and 


irreps for ^^^Tb 


label 


{\,fJ,)j,ST 


(A^, Unj^TT 


(Ajy, fJ'uj^U 


t 


(30, 4) i 


(10,4)1 


(20,0)0 





(28,8)1 


(10,4)1 


(18,4)0 


•k 


(30,4)1 


(11,2)1 


(18,4)0 


< 


(31,2)1 


(11,2)1 


(20,0)0 


• 


(31,2)1 


(10,4)1 


(18,4)0 


* 


(29,6)1 


(ll>2)i 


(18,4)0 





(29,6)1 


(10,4)1 


(18,4)0 


X 


(30,4)1 


(10,4)1 


(18,4)0 


> 


(31,2)1 


(11,2)1 


(18,4)0 



Table 4 
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Figure Captions 



Fig. 1: Experimental (left hand side) and theoretical (right hand side) ener- 
gies for the first four bands in ^^^Tb. 

Fig. 2: Experimental (left hand side) and theoretical (right hand side) ener- 
gies calculated without single-particle energies. 

Fig. 3: Energy bands in ^^^Tb. The first and tenth columns depict the 
experimental spectra. The second column shows the spectra obtained using 
Hamiltonian (0). The next seven columns show the same spectra obtained 
using Hamiltonian (|^) but with one term neglected, which is identified below 
each column: the third column does not include pairing, the fourth column 
does not include proton pairing, the fifth column does not include neutron 
pairing, the sixth column does not include the Kj term, the seventh column 
does not include the J^ term, the eighth eighth column does not include the 
asymmetry term, and finally, the ninth column does not include any of the 
three rotor terms. Insert (a) shows the 3/2^ band, insert (b) the 5/2^ band 
and insert (c) the 1/2^ band. 

Fig. 4: ^^^Tb theoretical energies and B{E2; J —>■ {J — 1)} transition 
strengths (in units e^6^ x 10~^). The latter are written close to the arrows 
indicating each transition. 

Fig. 5: ^^^Tb theoretical energies and B{E2] J —>■ {J — 2)} transition 
strengths, with the same convention as for Fig. 4. 

Fig. 6: SU(3) components of the calculated eigenstates in the first four 
bands of ^^^Tb. The angular momentum is listed in the horizontal axis. The 
symbols refer to the first column in Table 4. Insert (a) shows the 3/2^^ band, 
insert (b) the 5/2^ band, insert (c) the 1/2^^ band and insert (d) the 3/22 
band. 
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